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• Abstract 



In this report, the explicit probability density functions of the random Euclidean distances associated 
with equilateral triangles are given, when the two endpoints of a link are randomly distributed in 1) the 
same triangle, 2) two adjacent triangles sharing a side, 3) two parallel triangles sharing a vertex, and 4) 
two diagonal triangles sharing a vertex, respectively. The density function for 1) is based on the most 
recent work by Basel [1]. 2)— 4) are based on 1) and our previous work in 0, Q. Simulation results 
show the accuracy of the obtained closed-form distance distribution functions, which are important in 
ly-j . the theory of geometrical probability. The first two statistical moments of the random distances and the 

polynomial fits of the density functions are also given in this report for practical uses. 
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Index Terms 
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I. The Problem 

Define a "unit triangle" as the equilateral triangle with side length 1. Picking two points 
uniformly at random from the interior of a unit triangle, or between two adjacent unit triangles 
sharing a side or a vertex, the goal is to obtain the probabilistic density function (PDF) of the 
random distances between these two endpoints. There are four different cases \ab\, \pq\, \ef\ and 
\gh\, depending on the geometric locations of these two random endpoints, as shown in Fig. [Q 
The next section gives the explicit PDFs for these cases. 
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Fig. 1: Random Points and Distances Associated with Equilateral Triangles. 



II. Distance Distributions Associated with Equilateral Triangles 
A. \ab\: Distance Distribution within an Equilateral Triangle 

The author of LI J obtained the chord length distribution function for any regular polygon. From 
this result, HI further derived the density function for the distance between two uniformly and 
independently distributed random points in the regular polygon. Although the methods used were 
elementary, this work can be considered as a major breakthrough in Geometrical Probability, 
which also helps us verify the distance distribution in a regular hexagon (3). 

Following are the notations used in [|T||: V n ,r is the regular polygon with n vertices and with 
a circumscribed circle of radius r; l k is the distance between vertices, given by 

/ fc = 2rsin — , (1) 

n 

where k = 0, 1, K and K = [^rpj • L denotes the perimeter and A the area of V n y. 

it 1 2ir 

L = 2nr sin — and A = —nr 2 sin — . (2) 
n 2 n 



Denote the chord length distribution derived in [Q]| as f(s) for P nr , and the density function 
for the distance between two random points in V n>T as gp^d), the relationship between these 
two functions is as follows according to H|: 



0) 



According to this relationship and the derived chord length distribution f(s), the density 
function of random distances in an equilateral triangle, g Dl (d), is a special case in flU when 

n = 3, r 
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g Dl (d) = 4d< 



(Ad 2 + 6) sin^ 1 |f + ^2 - ^g) d 2 + 6V4d2 - 3 
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The corresponding cumulative distribution function (CDF) is 




Goi(d) = 2 ^ 



4d 2 



< d < 



V3 



Y < d < 1 



v/3 
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otherwise 



(4) 
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(d 2 + 3) sin- f d + (*f + l) y^Tf + (l - d^ 
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< d < 1 

d > 1 
(5) 



5. \pq\: Distance Distribution between Two Adjacent Equilateral Triangles Sharing a Side 

Given the result above, and the result obtained by us in |[2l for the distance distribution in a 
rhombus, further results of the random distances between two equilateral triangles can be derived. 
In Fig. [Q rhombus OABC can be decomposed into two congruent, adjacent equilateral triangles 
AOAB and AOBC. Picking points uniformly at random from the interior of this rhombus, then 



the points are equally likely to fall inside any one of these two triangles. Therefore, given the 
location of one endpoint of a random link, the second endpoint falls inside the same triangle as 
the first one with probability | (such as \ab\), and with probability | falls inside the adjacent 
triangle (such as \pq\). 

Suppose rhombus OABC in Fig. \T\ has a side length of 1, then the distribution of \ab\ is 
known as gD x (d) in © above. Denote the distribution of \pq\ as gjj A (d). The probability density 
function of the random distances between two uniformly distributed points that are both inside 
the same rhombus is /^(d) (see (1) m [El)- From the reasoning in the previous paragraph, 
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(6) 



and we have 



9D A (d) = 2f Dl (d) -g Dl {d). 



(7) 



Therefore, the probability density function of the random distances between two uniformly 
distributed points, one in each of the two adjacent unit triangles that are sharing a side, is 



(l + if ) Sln -.^ + (i* _ I) rf. _ fv ^3 



vsv 1 -*)-" 1 ^ + (573 -§WVB^3 
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< d < 



< d < 1 



1 < d< y/3 
otherwise 



(8) 



The corresponding CDF is 



GdM = 2< 



+ f d 3 + f z d 2 
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-(*I + 1 )' i2 -J 
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d < 
< d < & 



^ < d < 1 



1 < d < y/3 



1 d > V3 

(9) 

Note that although unit triangles are assumed in dU)-©, the distance distribution functions 
can be easily scaled by a nonzero scalar, for equilateral triangles of arbitrary side length. For 
example, let the side length of such triangles be s > 0, then 



G sD (d) = P(sD <d) = P(D < 



G 



d, 



Therefore, 



g sD (d) = <? D (- 

s s 



1 d 

~9d{~) 



(10) 



C. \ef\: Distance Distribution between Two Parallel Equilateral Triangles Sharing a Vertex 

This case corresponds to the random distance \ef \ in Fig. [TJ Here four unit triangles AOAB, 
AOBC, AOCD and ABCG together create a larger equilateral triangle AAGD with side 
length 2. According to (flOl) , the density function of distance distribution inside triangle AAGD 
is 92Di(d) = |<7r>i(f)» as s = 2. On the other hand, if we look at the two random endpoints of 
a given link inside the large triangle, they will fall into one of the two following cases: i) one 
of the endpoints falls inside one of the three unit triangles on the boarder of the large triangle, 
such as AOAB, AOCD or ABCG, with probability |; ii) one of the endpoints falls inside the 
unit triangle AOBC in the middle, with probability |. Each of these two cases includes several 
more detailed sub-cases as follows: 



Case i) Given the location of the first endpoint, the second endpoint will fall inside the same 
triangle as the first one (such as \ab\) with probability \, fall inside the adjacent triangle 
sharing a side (such as \pq\) with probability \, and fall inside one of the parallel triangles 
sharing a vertex (such as |e/|) with probability \. 

Case ii) When the location of the first endpoint is in AOBC, the second endpoint will fall inside 
the same triangle with probability \, and fall inside one of the adjacent triangles sharing 
a side with probability |. 
Denote the density function of random distance \ef\ as g Dp (d), we have the following 
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4 WO + \9D A { d ) + \gD P ( d ) 
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( d ) + ±9dM 
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Hence, 



8 2 4 d 

9dM = 3^2D : (rf) - 3 WO - 9D A ( d ) = 3 Wt^ 



-g Dl (d)-9D A (d)- (12) 



Therefore, the probability density function of the random distances between two uniformly 
distributed points, one in each of the two parallel unit triangles that are sharing a vertex, is 



g Dp (d) = M< 



(f - l) sin" 1 *f + d 2 - V4rf 2 - 3 - |d + ^ + 1 
*(* + 2 )-- 1 ^+(|-^) rf2 + 4 v / ^3 



4 I d 

V3 



8 J _ _8ti_ 

3" 3^ 



< d< ^ 



f <^<1 



V3 
2 

1 < d < V3 



y/3<d<2 
otherwise 



(13) 



The corresponding CDF is 
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(14) 



D. \gh\: Distance Distribution between Two Diagonal Equilateral Triangles Sharing a Vertex 

This case corresponds to the random distance \gh\ in Fig. Q] Here a regular hexagon is divided 
into six unit triangles. Looking at the two random endpoints of a given link inside the hexagon, 
the first endpoint can fall inside any one of the six triangles, and the second endpoint will i) 
fall inside the same triangle as the first one (such as \ab\) with probability |; ii) fall inside 
the adjacent triangle sharing a side (such as \pq\) with probability ~; iii) fall inside the parallel 
triangle sharing a vertex (such as \ef\) with probability |; iv) fall inside the diagonal triangle 
sharing a vertex (such as \gh\) with probability I. 

The density function of the random distances within a regular hexagon has been derived in (3]|, 
and we denote it as fn(d). Also denote the density function of random distance \gh\ as go u {d), 
we have 

/ir(<0 = a9Di{d) + o9D A ( d ) + IdDpid) + \g Du {d), (15) 

6 6 

or. 

g DD (d) = 6f H {d) - \g Dl {d) + 2g DA (d) + 2g Dp (d)] . (16) 



Therefore, the probability density function of the random distances between two uniformly 
distributed points, one in each of the two diagonal unit triangles that are sharing a vertex, is 
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The corresponding CDF is 
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III. Verification by Simulation 

Figure [2] plots the probability density functions, as given in ©, ([8]), (fl"3l) and (flTl) of the four 
random distance cases shown in Fig. [IJ Figure [3] shows a comparison between the cumulative 



9 



2 
1.8 
1.6 
1.4 
1.2 

u_ 

O 1 
D_ 

0.8 
0.6 
0.4 
0.2 


0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
Distance - d 

Fig. 2: Distributions of Random Distances Associated with Equilateral Triangles. 
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Fig. 3: Distribution and Simulation Results for Random Distances Associated with Equilateral 
Triangles. 

distribution functions (CDFs) of the random distances, and the simulation results by generating 
1, 000 pairs of random points with the corresponding geometric locations. Figure [3] demonstrates 
that our distance distribution functions are very accurate when compared with the simulation 
results. 
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TABLE I: Moments and Variance — Numerical vs Simulation Results 



Endpoint Geometry 


nni./L'!.,. 
rDrvMm 


M D 


M D 


VarD 


Within a 
Single Triangle 


9D l {d) 


0.3647918433s 


0.1666666667s 


0.0335935777s^ 


Sim 


0.3636606517s 


0.1654245278s 


0.0331164521s^ 


Between Two 
Adjacent Triangles 


9D A (d) 


0.6599648287s 


0.5s 


0.0644464249s^ 


Sim 


0.6597703174s 


0.4991102154s 


0.0637684635s z 


Between Two 
Parallel Triangles 


g Dp (d) 


1.0423971067s 


1.1666666667s 


0.0800749386s 2 


Sim 


1.0423259678s 


1.1655336144s 


0.0791204824s^ 


Between Two 
Diagonal Triangles 


9dM 


1.1880379828s 


1.5s 


0.0885657513s z 


Sim 


1.1896683303s 


1.5048850531s 


0.0896042750s^ 



IV. Practical Results 

A. Statistical Moments of Random Distances 

The distance distribution functions given in Section [n] can conveniently lead to all the statistical 
moments of the random distances associated with equilateral triangles. Given gD x {d) in (HJ), for 
example, the first moment (mean) of d, i.e., the average distance within a unit triangle, is 

rl 13 

Mj£ = / xg Dl (x)dx = - + — ln(3) w 0.3647918433, 
1 Jo 5 20 



and the second raw moment is 



M Di = J x 2 g Dl (x)dx = -, 



from which the variance (the second central moment) can be derived as 

0.0335935777. 



Var Dl = M^- ' 



When the side length of the unit triangle is scaled by s, the corresponding first two statistical 
moments given above then become 

MjS = 0.3647918433s, = - and Var D , = 0.0335935777s 2 . (19) 

i i g 

Table U lists the first two moments and the variance of the random distances in all four cases 
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TABLE II: Coefficients of the Polynomial Fit and the Norm of Residuals (NR) 



PDF 


Polynomial Coefficients 


NR 


( A\ 


10 10 x [0.006410 - 0.062752 0.283869 - 0.787308 1.497828 
— z.lwiyyo z.lboozo — l.fzlwo4 l.UbooUo — U.fc>14obo U.lyoo4U 
-0.056449 0.012613 -0.002124 0.000263 - 0.000023 0.000001 

0] 


U.UUzo4b 


I A\ 

9dM 


10 8 x [-0.000192 0.003514 - 0.029768 0.154523 - 0.549699 

T /lOm /II O A HOI C7H A 7HQH77 /I OflOOQO O 01 AC\CP. 

1.559704 - 0.636497 0.194663 - 0.043503 0.006952 - 0.000722 
0.000047 - 0.000002 0] 


0.086864 


9d p {a) 


10 7 x [-0.000069 0.001559 - 0.016145 0.101938 - 0.439383 
l.o/UoJ3y — O.ZUZ4U1 O.rl4zlo — /.o/oy*/ o.410o4U — t>.t>y/oU/ 
4.441839 - 2.155120 0.781941 - 0.206851 0.038480 - 0.004776 
0.000365 - 0.000015 0] 


0.105657 


9D D (d) 


10 7 x [0.000023 -0.000530 0.005579 - 0.035274 0.150448 
-0.460164 1.045986 - 1.805145 2.394131 - 2.453508 1.942627 
-1.182240 0.547275 -0.189550 0.047945 - 0.008552 0.001022 
0.000076 - 0.000003 0] 


0.075017 



given in Section HO It also gives the corresponding simulation results for verification purposes. 



B. Polynomial Fits of Random Distances 

Table HH lists the coefficients of the degree-20 polynomial fits of the original PDFs given in 
Section HH from d 20 to d°, and the corresponding norm of residuals. Figure ISa)-(d) plot the 
polynomials listed in Table [IT] with the original PDFs. From the figure, it can be seen that all 
the polynomials match closely with the original PDFs. These high-order polynomials facilitate 
further manipulations of the distance distribution functions, with a high accuracy. 

V. Conclusions 

In this report, we gave the closed-form probability density functions of the random distances 
associated with equilateral triangles. The correctness of the obtained results has been validated 
by simulation. The first two statistical moments and the polynomial fits of the density functions 
are also given for practical uses. 
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(b) Between two Adjacent Triangles Sharing a Side 




(c) Between two Parallel Triangles Sharing a Vertex 




(d) Between two Diagonal Triangles Sharing a Vertex 



Fig. 4: Polynomial Fit of the Distance Distribution Functions Associated with Equilateral 
Triangles. 
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